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The nuclear level densities of '"* '"Sn and the y-ray strength functions of iifi iisawsi^ below the neutron 
separation energy are extracted with the Oslo method using the (^He, ay) and ('He,'He'7) reactions. The level 
density function of "'Sn displays step-like structures. The microcanonical entropies are deduced from the level 
densities, and the single neutron entropy of "'Sn is determined to be (1.7±0.2) kg. Results from a combinatorial 
model support the interpretation that some of the low-energy steps in the level density function are caused by 
neutron pair-breaking. An enhancement in all the y-ray strength functions of "^""'Sn, compared to standard 
models for radiative strength, is observed for the y-ray energy region of (4-11) MeV. These small resonances 
all have a centroid energy of 8.0(1) MeV and an integrated strength corresponding to 1.7(9)% of the classical 
Thomas-Reiche-Kuhn sum rule. The Sn resonances may be due to electric dipole neutron skin oscillations or to 
an enhancement of the giant magnetic dipole resonance. 

PACS numbers: 21.10.Ma, 24.10.Pa, 24.30.Gd, 27.60.+j 



I. INTRODUCTION 

The level density and the y-ray strength function are av- 
erage quantities describing atomic nuclei. They are impor- 
tant for many aspects of fundamental and applied nuclear 
physics, including calculations of nuclear properties, like re- 
action cross sections. Such cross sections are used for calcula- 
tions in, e.g., reactor physics and nuclear waste management, 
and of nuclear reaction rates in astrophysics for modeling of 
nucleosynthesis in stars. 

The nuclear level density of nuclei is defined as the number 
of levels per unit of excitation energy. The entropy and other 
thermodynamic properties may also be determined from the 
level density. Structures in the level density are expected to 
be due to shell gaps, breaking of nucleon Cooper pairs, and/or 
changes in the nuclear shape. 

In the majority of previous experiments, the level density is 
measured either only at low energy by direct counting (con- 
ventional spectroscopy), or at higher energy around the neu- 
tron/proton separation energies (nuclear resonance measure- 
ments). 

The y-ray strength function may be defined as the reduced 
average transition probability as a function of y-ray energy. 
This quantity characterizes average electromagnetic proper- 
ties of excited nuclei. The strength function reveals essential 
information about the nuclear structure. Electric transitions 
are mostly influenced by the proton charge distribution, while 
magnetic transitions are also affected by the neutron distribu- 
tion due to the magnetic dipole moment of the neutron. The 
shape and softness of the nuclear surface are other important 
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factors for the nuclear response to electromagnetic radiation. 

The large number of stable isotopes in Sn makes the ele- 
ment suitable for systematic studies. This paper presents the 
level densities of "^•"''Sn and the y-ray strength functions 
of ii6.iis,ii9(5jj fQj. energies in the quasi-continuum below the 
neutron separation energy. The measurements have been per- 
formed at the Oslo Cyclotron Laboratory (OCL). A combina- 
torial model is also used in order to study, e.g., the origin of 
the level density steps and the impact of collective effects. 

The "'^■"''Sn results are compared with earlier OCL studies 
on other isotopes. In Ref 1 1 1, the level density functions of 
ii6,ii7|gjj were shown to display steps that are much more dis- 
tinctive than previously measured for other mass regions. The 
steps were interpreted as neutron pair-breaking. In Ref. [21, 
a resonance-like structure in the y-ray strength function was 
measured below the neutron threshold in "^Sn. 

The experimental set-up and the data analysis are briefly 
described in Sec. II. The normalized experimental results for 
level density and entropy are presented in Sec. III. Section IV 
discusses the nuclear properties extracted from the level den- 
sity with the combinatorial model. Section V presents the nor- 
malized experimental y-ray strength functions. Conclusions 
are drawn in Sec. VI. 



II. EXPERIMENTAL SET-UP AND DATA ANALYSIS 

The self-supporting "'Sn target was enriched to 93.2% 
and had a mass thickness of 1.6 mg/cm^. For three days 
the target was exposed to a 38-MeV ■'He beam with an av- 
erage cuiTent of ~ 1 .5 nA. The reaction channels studied were 
"^Sn(^He,3He'y)"'^Sn and "'^Sn(3He, Q-y)"'*Sn. 

Particle-y coincidences were recorded with 7 collimated Si 
particle A/i - E telescopes and 26 coIUmated Nal(TI) y-ray 
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detectors. The AE and E detector thicknesses were about 
140 yum and 1500 fim, respectively. These detectors were 
placed at 45° with respect to the beam axis. The Nal detectors 
are distributed on a sphere and constitute the CACTUS multi- 
detector system |3 1. The total solid-angle coverages out of An 
were approximately 1.3% for the particle detectors and 16% 
for the y-ray detectors. 

In the data analysis, the measured ejectile's energy is trans- 
formed into excitation energy of the residual nucleus using 
reaction kinematics. The y-ray spectra for various initial ex- 
citation energies are unfolded with the known response func- 
tions of CACTUS and the Compton subtraction method |3|. 
The Compton subtraction method preserves the fluctuations 
in the original spectra without introducing further, spurious 
fluctuations. 

The first generation y-ray spectra are extracted from the un- 
folded total y-ray spectra, by the subtraction procedure de- 
scribed in Ref. lIH. The main assumption is that the y-decay 
from any excitation energy bin is independent of the method 
of formation - whether it is directly formed by a nuclear reac- 
tion or indirectly by y-decay from higher lying states follow- 
ing the initial reaction. 

The first generation y-ray spectra are arranged in a 2- 
dimensional matrix P{E, Ey). The entries of P are the proba- 
bilities P{E, Ey) that a y-ray of energy Ey is emitted from an 
energy bin of excitation energy E. This matrix is used for the 
simultaneous extraction of the y-ray strength function and the 
level density function. 

The first generation matrix P is factored into the level den- 
sity function p and the radiative transmission coefficient T 
L5J: 



P(E,Ey)ccT(Ey)p(E-Ey). 



(1) 



The factorization of P into two components is justified for 
nuclear reactions leading to a compound state prior to a sub- 
sequent y-decay |6|. Equation (jTJ may also be regarded as a 
variant of Fermi's golden rule: The decay rate is proportional 
to the density of the final state and the square of the matrix 
element between the initial and final state. The factorization 
is performed by an iterative procedure where the independent 
functions p and T are adjusted until a global minimum with 
the experimental P(E, Ey) is reached. 

As shown in Eq. ([T]l, the transmission coefficient is assumed 
to be a function of only Ey, in accordance with the general- 
ized form of the Brink- Axel hypothesis |7 8 1. This hypothesis 
states that a giant electric dipole resonance, and all other col- 
lective excitation modes, may be built on any excited state and 
still have the same properties as the one built on the ground 
state. Hence, the transmission coefficient is independent of 
excitation energy. 

Equation ([T]l determines only the functional forms of p and 
T. The entries of P are invariant under the following transfor- 
mations ||5|: 

piE-Ey) = Aexp[a[E-Ey)]p{E-Ey), (2) 
t(Ey) = Bexp[aEy)T(Ey). (3) 

The final step of the Oslo method is to determine the normal- 
ization parameters. The parameters A and B will define the 



absolute values of p and 7~, respectively, while a will define 
their common slope. 



III. LEVEL DENSITIES 

A. Normalization and experimental results 

The constants A and a in Eq. Q, which are needed to nor- 
malize the experimental level density p, are determined using 
literature values of the known discrete energy levels at low 
energy and of the level spacing D at the neutron separation 
energy S„, obtained from neutron resonance experiments. 

The normaUzation value p(S„) is calculated either from the 
5-wave level spacing D(){S„) or from the /?-wave level spac- 
ing Di(S„). The level spacings are taken from Refs. llQl lTOl . 
To establish an expression for the value of p(5„), it is neces- 
sary to assume models for the spin distribution g(E, I) and the 
spin cutoff parameter cr. We choose the back-shifted Fermi 
gas (BSFG) model with the original parameterization of von 
Egidy et al. Ill], because this parameterization gives the most 
appropriate normalization of these nuclei when comparing to 
other experimental measurements (see also Ref. [jj). Here, 
these functions are kept as the original Gilbert and Cameron 
expressions 1 14|, but with a redefined parameterization of the 
nucleus' intrinsic excitation energy U and the level density 
parameter a. The spin distribution is expressed as IfTTI : 



g(EJ) 



21+ 1 

2cr2 



exp[-(/+ 1/2)2 /2cr2], 



where / is the spin, and where the spin cutoff" parameter a-(E) 
is given by: 



cr^ = 0.0888 A^^'^ar 



(5) 



where A is the mass number of the isotope, and T is the nu- 
clear temperature given by T = yfUja. Here, the level den- 
sity parameter is defined as a - 0.21 A'^^^ Me V"', while the 
shifted excitation energy U is defined as U - E - £paii - Ci. 
The back-shift parameter is defined as Ci = -6.6 A^ ^^^ MeV. 
The pairing energy Ep^i^ is calculated from the proton and neu- 
tron pair-gap parameters: iipaii- = Ap-i-An. The pair-gap param- 
eters are evaluated from the even-odd mass differences found 
in Ref. [15 | according to the method of Ref. liT6l . 

Assuming this spin distribution and equal numbers of levels 
with positive and negative parity, the level density at 5„ may 
be expressed as, for i-wave neutron resonances iBI fTTl : 



2cr2 , 

Po(S„) = t;- ■ + l)exp 



-dr + ir 



2(r2 



+ I, exp 



2(r2 



(6) 
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and for p-wave resonances ifTTl : 



-dr-ir 



2cr2 



+ /, exp 



2cr2 



+ (/, + 1) exp 
+ (I, + 2) exp 



■Ut + lf 



2cr2 



20-2 



(7) 



where the spin cutoff parameter is evaluated at 5„, and where 
/, is the spin of the target. 

A higher p{S„) is obtained from the level spacing of Dq 
than of Di, according to calculations on both isotopes. As the 
highest value of the level density is presumed to be the best 
estimate. Do is chosen in the following. The input parameters 
and the resulting values of the normalization data piS„) are 
given in Tab. [l] 

The experimental data for the level densities are not ob- 
tained up to the excitation energy of S„. There is a gap, and 
the level density in the gap and below is estimated according 
to the level density prediction of the BSFG model with the pa- 
rameterization of von Egidy et al. ifTTI . This is a consistency 
choice in order to keep the spin distribution and the spin cutoff 
parameter the same as the ones used during the calculation of 
p{S „) based on the neutron resonance data. The BSFG level 
density, for all spins and as a function of excitation energy, is 
given by 



p{E) 



BSFG 



exp(2 Vaf7) 



A scaling parameter rf is applied to the BSFG formula, 



p{E) 



BSFG 



TJpiE) 



BSFG : 



(8) 
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in order to make its absolute value at 5 „ agree with the nor- 
malization value p(S„). We then get a level density interpo- 
lation that overlaps with the measurements, and to which the 
measurements are normalized. The values of rj are shown in 
Tab.lH 

Figure [l] shows the normalized level densities in '"*'"^Sn. 
The arrows indicate the regions used for normalization. As 
expected, the level densities of "^Sn and "^Sn are very simi- 
lar to those of ''^Sn and "^Sn 1 1 1, respectively. 

The figure also shows that the known discrete levels fW\ 
seem to be complete up to ~ 2 MeV in "^Sn and up to 
~ 3 MeV in "^Sn. Hence, our experiment has filled a re- 
gion of unknown level density from the discrete region and 
to the gap, approximately at (S„ - 1) MeV. Unlike "^Sn, the 
ground-state of the even-even nucleus "^Sn has no unpaired 
neutron, and accordingly it has fewer available states than 
"^Sn. Therefore, measuring all levels to higher excitation en- 
ergies by conventional methods is easier in "**Sn. 

An alternative interpolation method to describe the gap be- 
tween our measured data and the neutron resonance data based 
piS „) is the constant temperature (CT) model Ill4il . This ap- 
proximation gives 



P(E) 



1 




exp [(E - Eo) IT] 



(10) 
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FIG. 1: Normalized level densities of "'Sn (upper panel) and "*Sn 
(lower panel) as a function of excitation energy. Our experimental 
data are marked with filled squares. The dashed lines are the BSFG 
predictions that are used for interpolation, scaled to coincide with 
p(5„) (open squares), which are calculated from neutron resonance 
data. The solid lines represent the discrete level densities obtained 
from counting the known levels. The arrows indicate the regions 
used to normalize the absolute values and the slope. The energy bins 
are 360 and 240 keV/ch for '"^ '''Sn, respectively. 



where the "temperature" T and the energy shift £0 are treated 
as free parameters. Figure |2] shows a comparison of the 
CT model and the BSFG model as interpolation methods for 
"'^Sn. The small difference in the region of interpolation is 
negligible for the normalization procedure. 



B. Step-like structures 

In Fig. [T[ the level density of "^Sn shows a step-like 
structure superimposed on the general level density, which is 
smoothly increasing as a function of excitation energy. A step 
is characterized by an abrupt increase of level density within 
a small energy interval. The phenomenon of steps was also 
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TABLE I: Input parameters and the resulting values for the calculation of the normalization value p(S„), and the input parameters for the BSFG 
interpolation and the required values of the scaling parameter rj. 



Nucleus 


S„ Dq{S„) a Ci An Ap o-(5„) p{S„) 
(MeV) (eV) (MeV-') (MeV) (MeV) (MeV) (IQ-^MeV-') 


1 


"«Sn 
i"*Sn 


6.485 700(150) 13.43 -1.43 1.02 4.55 6.05(175) 
9.326 61(7) 13.33 -1.43 1.19 1.24 4.74 38.4(86) 


0.44 
0.59 




6 7 8 

Excitation energy E (IVIeV) 

FIG. 2: (Color online.) Comparison of the BSFG model (dashed 
line) and the CT model (solid line) as interpolation means for the 
level density of "^Sn. The arrows indicate the region of normal- 
ization. The parameters in the CT model (T = 0.86 MeV and 
£0 = - 1 .7 MeV) have been found from a least fit to the data points 
in this region and from matching p(S„), and the parameterization is 
not intended to be appropriate elsewhere. 



seen in '"'■"^Sn Q. 

Distinctive steps in ''^Sn are seen below ~ 4 MeV. They 
are, together with the steps of "^''^''Sn |T|, the most distinc- 
tive steps measured so far at OCL. This may be explained by 
Sn having a magic number of protons, Z - 50. As long as the 
excitation energy is less than the energy of the proton shell 
gap, only neutron pairs are broken. The steps are distinctive 
since no proton pair-breaking smears the level density func- 
tion. 

The steps are less pronounced for "^Sn than for "^Sn. 
This is in contradiction to what is expected, as "**Sn is an 
even-even nucleus without the unpaired neutron reducing the 
clearness of the steps in "^Sn. The explanation probably lies 
in poorer statistics for the (^He, a) reaction channel than for 
(^He,^ He'). To reduce the error bars, a larger energy bin is 
chosen for "^Sn, leading to smearing the data and less clear 
structures. 

Two steps in "^Sn are particularly distinctive: one at 
~ 1 .0 MeV and another at ~ 2.0 MeV, leading to bumps in the 
region around ~ (1.2-1.4) MeV and around ~ (2.2-2.6) MeV, 
respectively. The steps in ' '''Sn are found at approximately the 
same locations as in "^Sn jH. 



Also for ''^Sn, two steps were clearly seen for low excita- 
tion energy [IJ. The first of these is probably connected to the 
isotope's first excited state, at 1.29 MeV ifTSl . A similar step 
in"«Sn would probably also had been found connected to the 
first excited state, at 1.23 MeV 1 18 1, if the measured data had 
had better statistics. 

Microscopic calculations based on the seniority model in- 
dicate that step structures in level density functions may be 
explained by the consecutive breaking of nucleon Cooper 
pairs 1,191 . The steps for "^Sn in Fig. [T] are probable candi- 
dates for the neutron pair-breaking process. The neutron pair- 
breaking energy of "''Sn is estimated' to be 2An - 2.5 MeV, 
which supports neutron pair-breaking as the origin of the pro- 
nounced bump around ~ (2.2 - 2.6) MeV. 

However, if the applied values of the neutron pair-gap pa- 
rameters are accurate, the pronounced step at ~ 1.0 MeV in 
"^Sn and other steps below this energy are probably not due 
to pure neutron pair-breaking. They might be due to more 
complex structures, involving collective effects such as vibra- 
tions and/or rotations. In Sec. |IV C] the pair-breaking in our 
isotopes will be investigated further. 



C. Entropy 

In many fields of natural science, the entropy is used to 
reveal the degree of order/disorder of a system. In nuclear 
physics, the entropy may describe the number of ways the 
nucleus can aiTange for a certain excitation energy. Vaiious 
thermodynamic quantities may be deduced from the entropy, 
e.g., temperature, heat capacity and chemical potential. The 
study of nuclear entropy also exhibits the amount of entropy 
gained from the breaking of Cooper pairs. We would like to 
study the entropy difference between odd-A and even-even Sn 
isotopes. 

The microcanonical entropy is defined as 



5,(£) = /tBlnQ.(£), 



(11) 



where is the Boltzmann constant, which is set to unity to 
make the entropy dimensionless, and where ils{E) is the state 
density (multiplicity of accessible states). The state density is 



' The values of the neutron pair-breaking 2An and the proton pair-breaking 
2Ap for '"* '"Sn are estimated from the A^/p values in Tab. I except for 
An of '"Sn. We estimate the energy for breaking a neutron pair in '"Sn 
as the mean value of 2An of the neighbouring even-even nuclei, redefining 
its value to be 2A„ = 2.5 MeV. 
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proportional to the experimental level density p(E) by 



a(£)(xp(£). [2 </(£)) + !] , 



(12) 



where {1(E)} is the average spin within an energy bin of exci- 
tation energy E. The factor 2 {I{E)) + 1 is the spin degeneracy 
of magnetic substates. 

The spin distribution is not well-known, so we assume the 
spin degeneracy factor to be constant and omit it. Omitting 
this factor is firstly grounded by the spin being averaged over 
each energy bin, leading to only the absolute value of the state 
density at high excitation energies being altered, and not the 
structure. Secondly, the average spin {1(E)} is expected to be 
a slowly varying function of energy (see Sec. |IV| i. Hence, a 
"pseudo" entropy S i may be defined based only on the level 
density p(E): 



Si(E)^kBln 



P(E) 
Po 



(13) 



The constant po is chosen so that Si - in the ground 
state of the even-even nucleus "^Sn. This is satisfied for 
Po - 0.135 MeV"^ The same value of po is used for "^Sn. 

Figure [3] shows the experimental results for the pseudo en- 
tropies of '"^Sn. These pseudo entropy functions are very 
similar to those of "^'"^Sn 1 1], which is as expected from the 
general similarity of the level density functions of these iso- 
topes. 

We define the entropy difference as 



AS(E)^Si^-S,'^-^ 



(14) 



where the superscript denotes the mass number of the isotope. 
Assuming that entropy is an extensive quantity, the entropy 
diff'erence will be equal to the entropy of the valence neutron, 
i.e. the experimental single neutron entropy of '^Sn. 

For midshell nuclei in the rare-earth region, a semi- 
empirical study |l20l has shown that the average single nu- 
cleon entropy is AS - 1.7 ks- This is true for a wide range 
of excitation energies, e.g., both for 1 and 7 MeV. Hence for 
these nuclei, the entropy simply scales with the number of nu- 
cleons not coupled in Cooper pairs, and the entropy difference 
is merely a simple shift with origin from the pairing energy. 

Figurejsjalso shows the entropy difference AS of '"*'"^Sn, 
which are midshell in the neutrons only. Above ~ 3 MeV, the 
entropy difference may seem to approach a constant value. In 
the energy region where the entropy difference might be con- 
stant (shown as the dashed line in Fig.[3]l, we have calculated 
its mean value as A5 = (1.7 + 0.2) fcg. Within the uncertainty, 
this limit is in good agreement with the general conclusion of 
the above mentioned semi-empirical study ll20l . and with the 
findings for "^ "^Sn fl^. For lower excitation energy, how- 
ever. Fig. [3] shows that the entropy difference of '"^ "^Sn is 
not a constant, unlike the rare-earth midshell nuclei. Hence, 
the "^•"''Sn isotopes have an entropy difference that is more 
complicated than a simple excitation energy shift of the level 
density functions. 
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FIG. 3: Upper panel: Experimental pseudo entropies 5/ of "'Sn 
(filled squares) and "*Sn (open squares) as a function of excitation 
energy. Lower panel: The respective experimental entropy differ- 
ence, AS = S/"' - 5/"**, as a function of excitation energy. An 
average value of AS(E) = (1.7±0.2)fca is obtained from a;^f--fit 
(dashed line) to the experimental data above ~ 3 MeV. 



IV. NUCLEAR PROPERTIES EXTRACTED WITH A 
COMBINATORIAL BCS MODEL 

A simple microscopic model ||2TI - |231 has been developed 
for further investigation of the underlying nuclear structure 
resulting in the measured level density functions. The model 
distributes Bardeen-Cooper-Schrieffer (BCS) quasi-particles 
on single-particle orbitals to make all possible proton and neu- 
tron configurations for a given excitation energy E. On each 
configuration, collective energy terms from rotation and vi- 
bration are schematically added. Even though this is a very 
simple representation of the physical phenomena, this combi- 
natorial BCS model reproduces rather well the experimental 
level densities. As a consequence, the model is therefore as- 
sumed to be able to predict also other nuclear properties of the 
system. We are first and foremost interested in investigating 
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Quadrupole deformation 

FIG. 4: (Color online) The Nilsson level scheme, showing single- 
particle energies as functions of quadrupole deformation ei, for 
"*Sn, which has 62 = 0.111. The Nilsson parameters are set to 
K = 0.070 and fi = 0.48. The Fermi levels are illustrated as curled 
lines (A„ for protons, and Ay for neutrons). 



the level density steps, and in investigating the assumption of 
parity symmetry, used in the normalization processes of the 
Oslo method. 



A. The model 

The single-particle energies e^p are calculated from the 
Nilsson model for a nucleus with an axially deformed core 
of quadrupole deformation parameter £2- The values of the 
deformation parameters are 62 = 0.111 and 62 - 0.109 |24| 
for ^'^'"^Sn, respectively. Also needed for the calculation of 
the Nilsson energy scheme are the Nilsson parameters k and fi 
and the oscillator quantum energy between the main oscillator 
shells: hcjo = 41A"'^''. The adopted values are k - 0.070 and 
fu = 0.48 for both neutrons and protons and for both nuclei, 
in agreement with the suggestion of Ref. [25|. All input pa- 
rameters are listed in Tab.[ll] The resulting Nilsson scheme for 
'^**Sn is shown in Fig. [4] 

The parameter A represents the quasi-particle Fermi level. 
It is iteratively determined by reproducing the right numbers 
of neutrons and protons in the system. The resulting Fermi 
levels for our nuclei are listed in Tab. and illustrated for 
ii^Sn in Fig. 



The microscopic model uses the concept of BCS quasi- 
particles (26'\. Here, the single quasi-particle energies gqp are 
defined by the transformation of: 

*^qp = -^(esp -^f + A2. (15) 

The pair-gap parameter A is treated as a constant, as before, 
and with the same values. 

The proton and neutron quasi-particle orbitals are charac- 
terized by their spin projections on the symmetry axes Q.„ and 
Q.y, respectively. The energy due to quasi-particle excitations 
is given by the sum of the proton and neutron energies and of 
a residual interaction V: 

£qp(Q„Qv)= ^ eqp(ii;) + eqp(n:,) + y(Q;,Q;). (16) 

In the model, quasi-particles having Q's of different sign will 
have the same energy, i.e. one has a level degeneracy. Since no 
such degeneracy is expected, a Gaussian random distribution 
V is introduced to compensate for a residual interaction appar- 
ently not taken into account by the Hamiltonian of the model. 
The maximum allowed number of broken Cooper pairs in our 
system is 3, giving a total of 7 quasi-particles for the even-odd 
nucleus "^Sn. Technically, all configurations are found from 
systematic combinations. 

On each configuration, both a vibrational band and rota- 
tions are built. The energy of each level is found by adding 
the energy of the configuration and the vibrational and rota- 
tional terms: 

E = E^piQ.^, Q,) + nco^i^y + A,otR (R+l) ■ (17) 

The vibrational term is described by the oscillator quantum 
energy ^Wyib and the phonon quantum number v = 0, 1,2, ... 
The values of /jWyib are found from the 2^ and 3" vibrational 
states of the even-even nucleus and are shown in Tab. The 
last term of Eq. ( [T7] i represents the rotational energy. The 
quantity Arot - ^/2J" is the rotational parameter with J' 
being the moment of inertia, and R is the rotational quan- 
tum number. The rotational quantum number has the val- 
ues of R = 0, 1,2, 3... for the even-odd nucleus of "^Sn, and 
R -0, 2,4... for the even-even nucleus of '"^Sn. 

For low excitation energy, the value of the rotational param- 
eter Aiot is determined around the ground state Ags. At high 
energy, the rotation parameter is found from a rigid, elliptical 
body, which is ETll : 



^rioid — n • y^oj 

° 4M/;2(i +0.3162) 

Here, M is the mass and Ra the radius of the nucleus. For 
nuclei in the medium mass region, A ~ 50 - 70, the rotational 
parameter A^gid is obtained at the neutron separation energy, 
according to a theoretical prediction [28 1. We assume that 
Aiigid is obtained at the neutron separation energy also for our 
nuclei. The applied values of Ags and Angid are listed in Tab.[ll| 
The function Ay^t as a function of energy is estimated from a 
linear interpolation between these. 
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TABLE II: Input parameters used in the combinatorial BCS model, and the resulting values for the Fermi levels Ay (neutrons) and A„ (protons). 



Nucleus 


(MeV) (MeV) (MeV) (MeV) (MeV) (MeV) 


"«Sn 
'"^Sn 


0.109 0.070 0.48 8.34 0.200 0.0122 1.23; 2.32 48.9 44.1 
0.111 0.070 0.48 8.36 0.205 0.0124 1.23; 2.32 48.9 44.2 



B. Level density 



In Fig. |5] the level density functions calculated by our 
model are compared with the experimental ones. We see that 
the model gives a very good representation of the level densi- 
ties in the statistical area above 3 MeV for both isotopes. Not 
taking into account all collective bands known from literature, 
the model is not intended to reproduce the discrete level struc- 
ture below the pair-breaking energy. The model also succeeds 
in reproducing the bump around ~ (2.2 - 2.6) MeV for both 
isotopes, even though the onset of this bump in "^Sn appears 
to be slightly delayed. Above ~ 2.6 MeV, the step brings the 
level density to the same order of magnitude as the experi- 
mental values. 

According to Eq. ([8]l and the relation between the intrinsic 
excitation energy U and the pair-gap parameters Ap and A,, the 
log-scale slope of the level density function is dependent on 
the pair-gap parameters. Figure|5]shows that the model repro- 
duces well the slopes of the level densities for both isotopes. 
This supports the applied values of the pair-gap parameters. 



C. Pair-breaking 

The pair breaking process produces a strong increase in 
the level density. Typically, a single nucleon entropy of 
(1.6 - 1.7) represents a factor of ~ 5 more levels due to 
the valence neutron. Thus, the breaking of a Cooper pair rep- 
resents about 25 more levels. Pair-breaking is the most impor- 
tant mechanism for creating entropy in nuclei as function of 
excitation energy. 

The average number of broken Cooper pairs per energy bin, 
(NqpY is calculated as a function of excitation energy by the 
model, using the adopted pair-gap parameters as input values. 
All configurations obtained for each energy bin are traced, and 
their respective numbers of broken pairs are counted. The av- 
erage number of broken pairs is also calculated separately for 
proton and neutron pairs. The result for "**'"^Sn is shown in 

Fig- HI 

Figurep^shows that the first pair-breaking for both ' '^'"^Sn 
is at an excitation energy around (2.2 - 2.6) MeV. That energy 
corresponds to the pair-breaking energy plus the extra energy 
needed to form the new configuration. The figure also shows 
that, according to the model, the pair-breakings here are only 
due to neutrons. The step in the average number of broken 
pairs is abrupt, and this number increases from to almost 1. 
This means that there is a very high probability for the nucleus 
to undergo a neutron pair-breaking at this energy. Provided 
that our values of the pair-gap parameters are reasonable, so 
that this intense step in the number of broken neutron pairs in 




123456789 
Excitation energy E (MeV) 



10 



FIG. 5: Level densities of "'Sn (upper panel) and ' "*Sn (lower panel) 
as a function of excitation energy. The solid lines are the theoretical 
predictions of the combinatorial BCS model. The squares are our 
experimental data. 



the model corresponds to the distinctive step in level density 
at ~ 2.0 MeV in "^Sn (see Sec. IIIi, that step in level density 
is probably purely due to neutron pair-breaking. 

The increases in the average number of broken pairs are 
abrupt also for certain other excitation energies, namely 
around (5 - 6) MeV and (8 - 9) MeV, as shown in Fig. |5] 
Here, we predict increases of the numbers of levels, caused 
by pair-breaking, even though not necessarily visible with the 
applied experimental resolution. In between the abrupt pair- 
breakings, the number of broken pairs is almost constant and 
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FIG. 6: (Color online) The average number of broken quasi-particle 
pairs (a^,;,^ (solid line) as a function of excitation enegy for '"Sn 
(upper panel) and "*Sn (lower panel), according to the combinato- 
rial BCS model. Also shown is how this quantity breaks down into 
neutron pairs (dashed line) and proton pairs (dashed-dotted line). 



close to integers. Saturation has been reached, and signifi- 
cantly more energy is needed for the next pair-breaking. 

Neutron pair-breaking dominates over proton pair-breaking 
for the energies studied. Even though there is a large shell 
gap for the protons, breaking of proton pairs also occurs, but 
then only for energies above the proton pair-breaking energy 
of 2Ap plus the shell gap energy. According to Fig. |6j pro- 
ton pair-breaking contributes for excitation energies above 
3.5 MeV in both isotopes. An increased number of broken 
proton pairs at higher energies is expected to lead to the level 
density steps at high excitation energy being smeared out and 
becoming less distinctive, in accordance with the experimen- 



tal findings of Sec. Ill B 



Two effects due to the Pauli principle are notable in Fig.|6] 
In "^Sn compared to "**Sn, 1) the increases of the total aver- 
age number of broken pairs occur at higher energies; and 2) 
the average number of broken proton pairs is generally higher 



The explanation probably is that the valence neutron in "'Sn 
to some extent hinders the neutron pair-breaking. The pres- 
ence of the valence neutron makes fewer states available for 
other neutrons, due to the Pauli principle. Therefore, in ^'^Sn 
compared to "^Sn, more energy is needed to break neutron 
pairs, and for a certain energy, proton pair-breaking is more 
probable. Of course, an increase in the number of broken pro- 
ton pairs leads to a corresponding decrease in the number of 
broken neutron pairs. 



D. Collective effects 

We have made use of the model to make a simple estimate 
of the relative impact on the level density of collective effects, 
i.e., rotations and vibrations, compared to that of the pair- 
breaking process. The enhancement factor of the collective 
effects is defined as 



FcoiiiE) 



P(E) 



Pnon-coll(^) 



(19) 



where Pnon-coii is the level density function excluding collec- 
tive effects. 

Figure [T] shows the calculated level density with and with- 
out collective contributions from vibrations and from rota- 
tional bands for "^Sn. The model prediction is assumed to 
be reasonably valid above ~ 3 MeV of excitation energy. Ac- 
cording to these simplistic calculations, the enhancement fac- 
tor of collective effects F sharply decreases at the energies of 
the steps in the average number of broken quasi-particle pairs 
(see Fig.|6]). For ' ''^Sn, we find that F decreases for excitation 
energies of approximately 2.5 and 6 MeV, where the average 
number of broken quasi-particle pairs increase from ~ to 
~ 1, and from ~ 1 to ~ 2, respectively. For the energies stud- 
ied, the maximum value of F is about 10, found at £ ^ 6 
MeV. For '"*Sn, the enhancement factor would be less than 
for ''^Sn, since this nucleus does not have unpaired valence 
neutrons. 

As a conclusion, the collective phenomena of vibrations 
and rotations seem to have a significantly smaller impact on 
the creation of new levels than the nucleon pair-breaking pro- 
cess, which has an enhancement factor of typically about 25 
for each broken pair 



E. Parity asymmetry 



The parity asymmetry function a is defined as 



P+ -P- 
P+ +P- 



(20) 



where p+ is the level density of positive parity states, and p_ 
is the level density of the negative parity states. The values of 
a range from -1 to +1. A system with o- = -1 is obtained 
for p+ = 0, implying that all states have a negative parity. A 
system with a = has equally many states with positive as 
negative parity and is obtained for p_ = p+. 
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FIG. 7: (Color online) The impact of collective effects on the level 
density of "'Sn according to the combinatorial BCS model. The up- 
per panel shows experimental level density (data points) compared 
with model calculations with collective effects (solid line) and with- 
out collective effects (dashed line). The lower panel shows the cor- 
responding enhancement factor of the collective effects, F^^]] (linear 
scale). 

The Nilsson scheme in Fig. [4] shows that the single-particle 
orbitals both above and below the neutron Fermi level are a 
mixture of positive and negative parities. In addition, each of 
these states may be the head of vibrational bands, for which 
the parity of the band may be opposite of that of the band 
head. 

The parity asymmetry functions of "** "^Sn are drawn in 
Fig. [8] For energies below the neutron pair-breaking energy 
approximately at 2.5 MeV, the even-odd isotope has a parity 
asymmetry function with large fluctuations between positive 
and negative values, while the even-even isotope has positive 
parities. This is as expected when vibrational bands of oppo- 
site parity are not introduced. (The zero parity of "^Sn for 
certain low-energy regions is explained by the non-existence 
of energy levels.) 

Above the pair-breaking energies, the asymmetry functions 



FIG. 8: The parity asymmetry function a, according to the combina- 
torial BCS model, shown as a function of excitation energy for '"Sn 
(upper panel) and "*Sn (lower panel). 

begin to approach zero for both isotopes. This is also as ex- 
pected, since we then have a group of valence nucleons that 
will randomly occupy orbitals of positive and negative parity 
and on average give an a close to zero. Above ~ 4 MeV, the 
parity distributions of "^'"^Sn are symmetric with p+ ^ p_. 
This is a gratifying property, since parity symmetry is an 
assumption in the Oslo method normalization procedure for 
both the level density and the y-ray strength function (see 
Sec.|V|. 

F. Spin distribution 

The combinatorial BCS model determines the total spin / 
for each level from the relation: 

I{I+1) = R(R + 1)+ n^ + Qy, (21) 

from which the spin distribution of the level density p may be 
estimated. 
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FIG. 9: Comparison of spin distributions in '"*Sn at different excitation energies. Tlie open squares are calculated from the combinatorial BCS 
model in Eq. \2i\ , while the solid lines are the predictions of Gilbert and Cameron in Eq. (j4]l. The spin distributions are normalized to 1. 
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The resulting spin distribution is compared with the theo- 
retical spin distribution of Gilbert and Cameron in Eq. Q, 
using the same parameterization of the spin cutoff parameter 
as Eq. Q. Figure [9] shows the comparison for four different 
excitation energies: 5, 6, 7 and 8 MeV. The agreement is gen- 
erally good. Hence, the spin calculation in Eq. (21 1, and the 
assumption of the rigid rotational parameter A^gid obtained at 
S„, are indicated to be reasonable assumptions. We also note 
from the figure that the average spin, {!(£)}, is only slowly in- 
creasing with excitation energy, justifying the pseudo entropy 
definition introduced in Eq. (fT3]l. 



V. GAMMA-RAY STRENGTH FUNCTIONS 

A. Normalization and experimental results 

The y-ray transmission coefficient T, which is deduced 
from the experimental data, is related to the y-ray strength 
function / by 



XL 



(22) 



where X denotes the electromagnetic character and L the 
multipolariy of the y-ray. The transmission coefficient T is 
normalized in slope a and in absolute value B according to 
Eq. (pll. The slope was determined in Sec. Ill in the case of 

118,1 ftc^ 



^Sn, and in Ref. [1| in the case of "^Sn. The absolute 
value normalization is yet to be determined. This is done us- 
ing the literature values of the average total radiative width at 
the neutron separation energy, (ry(S„)Y which are measured 
for neutron capture reactions (n, y). 

The y-ray transmission coefficient T is related to the aver- 
age total radiative widthjTyiE, I, tt)^ of levels with energy E, 
spin / and parity k by {29] : 



(r,(£,/,7r)) = 



1 



2n p(E, I, jf) 



XL If,nf ^'^y 



dEyTxL(Ey)p(E - Ey, If, 71 f) . (23) 



The summations and integration are over all final levels of spin 
// and parity jif that are accessible through a y-ray transition 
categorized by the energy Ey, electromagnetic character X and 
multipolarity L. 

For i-wave neutron resonances and assuming a major con- 
tribution from dipole radiation and parity symmetry for all ex- 
citation energies, the general expression in Eq. ( [23] l will atS„ 
reduce to 

{Yy(S„,I,±l/2,7r,)) 

B r^" 

" ~A ,c J n — T dEyT{Ey)p{S„-Ey) 

I 

xY^g{S„-EyJ,±\l2 + J). (24) 



Here, /, and n, are the spin and parity of the target nucleus 
in the (n, y) reaction. Indeed, the results from the combina- 
torial BCS model in Sec. |IV| supports the symmetry assump- 
tion of the parity distribution. The normalization constant B 
in Eq. (j24| is determined |30| by replacing T with the ex- 
perimental transmission coefficient, p with the experimental 
level density, g with the spin distribution given in Eq. Q, and 
{VyiSn)) with its literature value. 

The input parameters needed for determinin g the normal- 
ization constant B for ^'^'^'^Sn are shown in Tab.llllland taken 
from Ref. [TOl. For "^Sn, the level spacing Do{S„) is not 
available in the literature. Therefore, p(S„) was estimated 
from systematics for the normalization of a in Ref. H]. The 



value of Do in Tab. Ill is estimated from piS„). Note that there 
was an error in the spin cutoff' parameters cr{S„) of "^'"^Sn 
in Refs. |[ll|2l. The impact of this correction on the normal- 
ization of level densities and strength functions is very small. 
Moreover, updated values of Do(5'„) and (ry{S„)^ are now 
available for "^Sn ifTOl . All the ne w no rmaUzation parameters 

The value of (ry(S„)^ 



for 



116,117 



Sn are presented in Tab. 



IV 



of ''^Sn is taken from the indicated value in Ref. ifTOl . 



TABLE III: Input parameters for normalization of the y-ray trans- 
mission coefficient T for '"* '"Sn. 



Nucleus 


/, Do(S„) (r,(5„)) 
{n) (eV) (meV) 


'"^Sn 


700 45 
1/2 61 117 



TABLE IV: New normalization parameters for ' 



^Sn. 



Nucleus 


a-{S„) Do(S„) p(S„) (r,(5„)) 
(eV) (10* MeV-') (meV) 


1 


"'Sn 
""Sn 


4.58 450(50) 9.09(2.68) 52 
4.76 59 40(20) 120 


0,43 
0,45 



The resulting y-ray strength functions of "^""^Sn are 
shown in Fig, [TO] For all isotopes, it is clear that there is a 
change of the log-scale derivate at Ey ~ 4 MeV, leading to 
a sudden increase of strength. Such an increase may indicate 
the onset of a resonance. The comparison in Fig, [TOjof the 
new "^Sn strength function with the earlier published one |i2l 
confirms that correcting the o-(S„) and the Do(S„) had only a 
minor impact on the normalization. 



j=-i 



Figure 1 1 shows the four normalized strength functions of 
ii6-ii9jgjj together They are all approximately equal except 
for '"^Sn, which has a lower absolute normalization than the 
others. This is surprising considering the quadrupole defor- 
mation parameter of ' '^Sn (e2 =0,111) being almost identical 
to that of '"'Sn (ei = 0,112) f24l. In the following, we there- 
fore multiply the strength of "**Sn with a factor of 1,8 to get it 
on the same footing as the others. The values of p(S„) and the 
scaling parameter r] (se e also Sec, |lII A"] ) of these Sn isotopes 
are collected in Tab,M For "'^Sn, the p(S „) may be expected 
to be larger, while 77 may be expected to be smaller It would 



12 



1 10-^ 



c 
_o 

u 
c 

3 



119, 



Sn 



t 



•ft* * 



I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ 




■ Present 
□ Published in [2] 



I I I I I I ^1 I I I I I I I I I I 



■ I ■ ■ ■ ■ I ■ ■ ■ ■ i ■ ■ ■ ■ I ■ ■ ■ ■ I 



1 10-^ 



118, 



Sn 



c 
_o 

u 
c 

3 



=io-« 

<u 



It 



I I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I I 



116, 



Sn 



J 



I ■ ■ ■ ■ I ■ ■ ■ ■ 



I I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I 



■'0 01234567891 23456789 
Y-ray energy E (MeV) y-ray energy E (MeV) 



FIG. 10: Normalized y-ray strength functions as functions of y-ray energy. The upper panels show '"Sn (left) and both present and previously 
published versions of "'Sn (right). The energy bins 120 and 240 keV/ch for "' "'Sn, respectively. The lower panels show "*Sn (left) and 
"*Sn (right) with the energy bins of 360 and 120 keV/ch, respectively. 



be desirable to remeasure both Do(S„) and (Vy^ for this iso- 
tope, since the apparent wrong normalization of the strength 
function of "^Sn depends on these parameters. 



TABLE V: TheFermigas approximation forp(5„)BSFG> the calculated 
p(5„), and the resulting scaling parameter rj, for "*""'Sn. 



Nucleus 


p(S ;,)bsfg 


P(5„) 


1 




(10" MeV- 


) (10" MeV-') 




'"Sn 


14 


6.05(175) 


0.44 


"«Sn 


65 


38.4(86) 


0.59 


'"Sn 


22 


9.09(2.68) 


0.43 


"^Sn 


89 


40.0(20.0) 


0.45 



B. Pygmy resonance 

Comparing our measurements with other experimental data 
makes potential resonances easier to localize. Experimental 



cross section data o-(Ey) are converted to y-ray strength f{Ey) 
through the relation: 



f(Ey) 



1 



cr(Ey) 



(25) 



Figure [T2| shows the comparison of the Oslo strength func- 
tions of '"'-"'^Sn with those of the photoneutron cross sec- 
tion reactions "*'"^Sn(y, n) from Utsunomiya et al. |13| and 
'''Sn('y, «) from Varlamov et al. 1311 . photoabsorption re- 



116-119 



Sn(7, x) from Fultz et al. 



116-118 



Sn(y, X) 



actions 

from Varlamov et al. f33l, and ''^-''**Sn(7, x) from Lepretre 
et al. [34] , Clearly, the measurements on "^•"'^Sn both from 
Oslo and from Utsunomiya et al. ifTSl independently indicate 
a resonance, from the changes of slopes. For "^'''*'Sn, the 
Oslo data clearly shows the presence of resonances. Hence, 
the resonance earlier observed in '^^Sn JJ) is confirmed also in 
Ii6,ii8,ii9igjj -pjjjg resonance will be referred to as the pygmy 
resonance. 

In order to investigate the experimental strength functions 
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FIG. 1 1 : (Color online) The four normalized strength functions of 
ii6-ii9gj^ shown together. 



farther, and in particular the pygmies, we have applied com- 
monly used models for the Giant Electric Dipole Resonance 
(GEDR) and for the magnetic spin-flip resonance, also known 
as the Giant Magnetic Dipole Resonance (GMDR). 

For the GEDR resonance, the Generalized Lorentzian 
(GLO) model 1351 is used. The GLO model is known to agree 
well both for low y-ray energies, where we measure, and for 
the GEDR centroid at about 16 MeV. The strength function 
approaching a non-zero value for low Ey is not a property spe- 
cific for the Sn isotopes, but has been the case for all nuclei 
studied at the OCL so far 

In the GLO model, the £1 strength function is given by 
03: 



1 



£1 



+ 0.7 



(£^^-£■£1^) + Ey^ (rKMF(Ey,Tf)j 

^KMpiEy - 0, Tf) 



(26) 



in units of MeV. 

Usually, Tf is interpreted as the nuclear temperature 
of the final state, with the commonly applied expression 
Tf - yfUja. On the other hand, we are assuming a con- 
stant temperature, i.e., the y-ray strength function is indepen- 
dent of excitation energy. This approach is adopted for con- 
sistency with the Brink- Axel hypothesis (see Sec. |ll]), where 
the strength function was assumed to be temperature indepen- 
dent. Moreover, we treat Tf as a free parameter in order to fit 
in the best possible way the theoretical strength prediction to 
the low energy measurements. The applied values of Tf are 
listed in Tab. ED 

TABLE VI: Parameters used for the theoretical y-ray strength func- 
tions of "*""'Sn. The value of Tf in '"*Sn has been found for the 
measured strength function multiplied by 1.8. 



Nucleus 


Eei 






Emi 


Lmi 


o'mi Tf 




(MeV) (MeV) 


(mb) 


(MeV) (MeV) 


(mb) (MeV) 


""Sn 


15.53 


4.81 


253.0 


8.34 


4.00 


0.963 0.40(1) 


"«Sn 


15.59 


4.77 


256.0 


8.36 


4.00 


0.956 0.40(1) 


'"Sn 


15.66 


5.02 


254.0 


8.38 


4.00 


1.04 0.46(1) 


""Sn 


15.68 


4.19 


266.0 


8.41 


4.00 


0.773 0.46(1) 



The Ml spin-flip resonance is modeled with the functional 
form of a Standard Lorentzian (SLO) l24l : 



fm°(Ey) 



1 



O" Ml^Ml''Ey 



') +Ey^rM\" 



(28) 



where the parameter Emi is the centroid energy, Fmi the width 
and (Tmi the cross section, of the GMDR. These Lorentzian 
parameters are predicted from the expressions in Ref. 1241 . 
with the results as shown in Tab. |Vl] 

In the absence of any established theoretical prediction 
about the pygmy resonance, we found that the pygmy is satis- 
factorily reproduced by a Gaussian distribution [2 J: 



fpyg(Ey) — C| 



1 



pyg 



V2^c 



exp 



{Ey ^pyg) 
2a- 2 

Z.U pyg 



(29) 



where Cpyg is the pygmy's normalization constant, Epy^ the 
energy centroid and cTpyg the standard deviation. These pa- 
rameters are treated as free. The total model prediction of the 
y-ray strength function is then given by: 



in units of MeV"^, where the Lorentzian parameters are the 
GEDR's centroid energy Ee\, width F^i and cross section 
(Tei ■ We use the experimental parameters of Fultz |32 1, shown 
in Tab.|Vl| The GLO model is temperature dependent from the 
incorporation of a temperature dependent width T^mf- This 
width is the term responsible for ensuring the non-vanishing 
GEDR at low excitation energy. It has been adopted from the 
Kadmenskii, Markushev and Furman (KMF) model I.36J and 
is given by: 



FKMFiEy, Tf) = ^ [Ey^ + 4n%^) 



(27) 



/to 



/eI + /mi + /p; 



(30) 



By adjusting the Gaussian pygmy parameters to make the 
best fit to the experimental d ata o f "^""^Sn, we obtained the 

The pygmy fit of "'Sn is 



Vll 



values as presented in Tab. 
updated and corresponds to the present normalization of the 
strength function. The pygmy fit also gave an excellent fit for 
''^Sn. For "**'"'^Sn, it was necessary to slightly reduce the 
values of Tf and cTpyg. The similarity of the sets of parame- 
ters for the different nuclei is gratifying. As is seen in Fig. [12] 
these theoretical predictions describe the measurements rather 
well. 
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FIG. 12: (Color online) Comparison of theoretical predictions including pygmy fits with experimental measurements for "*""'Sn. The total 
strengths (solid lines) are modeled as Gaussian pygmy additions to the GLO (El + Ml) baselines. The SLO (El + Ml) baselines are also 
shown, failing to reproduce the measurements for low Ey. The arrows indicate the neutron separation energies S„. 

Upper, left panel: Comparison of theoretical predictions of "'Sn with the Oslo measurements, '"Sn{y,n) from Utsunomiya et al. |13|, 
"'Sn(7, x) from Fultz ef a/. |32|,and "''Sn(7, «) from Varlamov et al. \2>\ \. 

Upper, right panel: Comparison of theoretical predictions of "'Sn with the Oslo measurements, "'Sn(7, «) from Utsunomiya et al. |13|, 
"■'Sn(7, x) from Fultz e/ a/. |32|, "'Sn(7, x) from Varlamov ef a/. |,33J , and '"Sn(7, x) from Lepretre e/ aZ. |34|. 

Lower, left panel: Comparison of theoretical predictions of "*Sn with the Oslo measurements multiplied with 1.8 (filled squares) (the mea- 
surements with the original normalization are also included as open squares), "^Sn(7, fi) from Utsunomiya et al. [[13J, "*'Sn(7, x) from Fultz 
et al. [32J, "**Sn(7, x) from Varlamov et al. |'33], and "*Sn(7, x) from Lepretre et al. 1341 . 

Lower, right panel: Comparison of theoretical predictions of "*Sn with the Oslo measurements, "*Sn(7,«) from Utsunomiya et al. |13|, 
"''Sn(7, x) from Fultz et al. 132J, "^Sn(7, x) from Varlamov et al. (331, and "*Sn(7, x) from Lepretre et al. l34l . 



The pygmy centroids of all the isotopes are estimated to be 
around 8.0(1) MeV. It is noted that an earlier experiment by 
Winhold et al. |37| using the (y, n) reactions determined the 
pygmy centroids for ''^'"^Sn to approximately 7.8 MeV, in 
agreement with our measurements. 

Extra strength has been added in the energy region of 
~ (4-11) MeV. The total integrated pygmy strengths are 
30(15) MeV-mb for all four isotopes. This constitutes 1.7(9)% 
of the classical Thomas-Reiche-Kuhn (TRK) sum rule, as- 
suming all pygmy strength being E\. Even though these res- 
onances are rather small compared to the GEDR, they may 



have a non-negligible impact on nucleosynthesis in super- 
novas 1381 . 

If one does not multiply the strength function of '"^Sn by 
1.8 for the footing equality, then the pygmy of "**Sn becomes 
very different from those of the other isotopes, and the total 
prediction is not able to follow as well the measurements for 
low Ey. A pygmy fit of the original normalization does how- 
ever give: Tf = 0.28(2) MeV, Cpyg = 1.8(6) ■ 10-^MeV-^ 
£pyg = 8.0(2) MeV and cr^y^ = 1.0(1) MeV. This rep- 
resents a smaller pygmy, giving an integrated strength of 
17(8) MeV-mb and a TRK value of 1.0(5)%. 
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TABLE VII: Empirical values of "*"'"Sn pygmies Gaussian param- 
eters, the integrated pygmy strengths and the TRK values. For "**Sn, 
the values have been found from fitting to the measured strength 
function multiplied by 1.8. 



Nucleus 


C F (T 
(lO-'MeV^-) (MeV) (MeV) 


Int. strength TRK value 
(MeV-mb) (%) 


"«Sn 
i"*Sn 
'"Sn 
'"^Sn 


3.2(3) 8.0(1) 1.2(1) 
3.2(3) 8.0(1) 1.2(1) 
3.2(3) 8.0(1) 1.4(1) 
3.2(3) 8.0(1) 1.4(1) 


30(15) 1.7(9) 
30(15) 1.7(9) 
30(15) 1.7(9) 
30(15) 1.7(9) 



The commonly applied Standard Lorentzian (SLO) was 
also tested as a model of the baseline and is included in 



Fig. 12 The SLO succeeds in reproducing the (y, x) data, but 
clearly fails for the low-energy strength measurements, both 
when it comes to absolute value and to shape. The same has 
been the case also for many other nuclei measured at the OCL. 
Therefore, we deem the SLO not to be adequate below the 
neutron threshold. 

Probably, the pygmies of all the Sn isotopes are caused by 
the same phenomenon. It is still indefinite whether the Sn 
pygmy is of £1 or Ml character. A clarification would be of 
utmost importance. 

Earlier studies indicate an £1 character of the Sn pygmy. 
Amongst these are the nuclear resonance fluorescence exper- 
iments (NRF) performed on "''■i^^Sn |39| and on ^^^'^^^Sn 
||40l, and the Coulomb dissociation experiments performed on 
i29-i32|5jj [41, 42|. If the Sn pygmy is of £1 character, it may 
be consistent with the so-called E\ neutron skin oscillation 
mode, discussed in Refs. |'43'-45 1. 

However, the possibility of an Ml character cannot be ruled 
out. Figure |4] shows that the Sn isotopes have their proton 
Fermi level located right in between the ^7/2 and ^9/2 orbitals, 
and their neutron Fermi level between /zii/2 and /19/2. Thus, 
an enhanced Ml resonance may be due to proton ^7/2 <-> ^9/2 
and neutron /111/2 <-» /19/2 magnetic spin-flip transitions. The 
existence of an Ml resonance in this energy region has been 
indicated in an earlier experimental study: proton inelastic 
scattering experiment on '^" '^"^Sn li46J . 



VI. CONCLUSIONS 

The level densities of "^ '^^Sn and the y-ray strength func- 
tions of ijaye been measured using the ("'He, ay) 



and (^He,^He'7) reactions and the Oslo method. 

The level density function of ' '''Sn shows pronounced steps 
for excitation energies below ~ 4 MeV. This may be explained 
by the fact that Sn has a closed proton shell, so that only neu- 
tron pairs are broken at low energy. Without any proton pair- 
breaking smearing out the level density function, the steps 
from neutron pair-breaking remain distinctive. The entropy 
has been deduced from the experimental level density func- 
tions, with a mean value of the single neutron entropy in ' '^Sn 
determined to (1 .7 ± 0.2) ks- These findings are in good agree- 
ment with those of '"" "^Sn. 

A combinatorial BCS model has been used to extract nu- 
clear properties from the experimental level density. The num- 
ber of broken proton and neutron pairs as a function of excita- 
tion energy is deduced, showing that neutron pair-breaking is 
the most dominant pair-breaking process for the entire energy 
region studied. The enhancement factor of collective effects 
on level density contributes a maximum factor of about 10, 
which is small compared to that of pair-breaking. The parity 
distributions are found to be symmetric above ~ 4 MeV of 
excitation energy. 

In all the ^'^ "^Sn strength functions, a significant en- 
hancement is observed in the energy region of Ziy ^ (4 - 1 1) 
MeV. The integrated strength of the resonances correspond to 
1 .7(9)% of the TRK sum rule. These findings are in agreement 
with the conclusions of earlier studies. 
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